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The dynamics of particles interacting by key-lock binding of attached biomolecules are studied 
theoretically. Examples of such systems include DNA-functionalized colloids as well as nanoparti- 
cles grafted with antibodies to cell membrane proteins. Depending on the coverage of the functional 
groups, we predict two distinct regimes separated by a percolation transition. In the localized regime 
at low coverage, the system exhibits a broad, power law like distribution of particle departure times. 
At higher coverage, there is an interplay between departure dynamics and particle diffusion. This 
interplay leads to a sharp increase of the departure times, a phenomenon qualitatively similar to 
aging in glassy systems. This diffusive regime is analogous to dispersive transport in disordered 
semiconductors: depending on the interaction parameters, the diffusion behavior ranges from stan- 
dard diffusion to anomalous, subdiffusive behavior. The connection to recent experiments and 
implications for future studies are discussed. 



Selective key-lock interactions are quintessential for 
biology. Over the past several years, they have 
also attracted substantial attention in the context of 
nanoscience. It is becoming common practice to at- 
tach biomolecules capable of key-lock binding to colloidal 
particles or other microscopic objects to achieve control- 
lable, specific interactions. Examples include nanopar- 
ticles functionalized with complementary single-stranded 
DNA (ssDNA) ([1[-[8|), or with antibodies to a particu- 
lar protein. The possible applications range from self- 
assembly of smart nanomaterials to biosensors and cell- 
specific drug delivery (|9(-[ll|). In this new class of sys- 
tems, the collective character of the binding may lead to 
non-trivial and often prohibitively slow dynamics. 

In this letter we report several remarkable results deal- 
ing with the dynamics of particles with reversible key- 
lock interactions. These results are of both conceptual 
and practical interest. In particular, we will demonstrate 
that depending on the coverage of the functional groups 
(e.g. ssDNA or proteins), the system exhibits two dis- 
tinct regimes separated by a percolation transition. At 
low coverage, there is a broad power law like distribution 
of departure times but no lateral diffusion. If the cover- 
age is sufficiently high, the overall particle dynamics is a 
result of the interplay between diffusion and desorption. 
The lateral motion is analogous to dispersive transport 
in disordered semiconductors: it may range from regu- 
lar diffusion with a renormalized diffusion coefficient, to 
anomalous, subdiffusive behavior. 

In the simplest version of our model, a single parti- 
cle interacts with a flat 2D surface via multiple key-lock 
binding (see Figure Q]) . At each position of the parti- 
cle, there are m key-lock bridges which may be closed 
or open, and there is a binding energy e for each of the 
key-lock pairs (the variation in e is neglected). Therefore, 
the m-bridge free energy plays the role of an effective lo- 
cal potential for the particle: U (m) = — fc^TmA, where 
A = log (1 + cxp [e/ksT]) [l2j. In a generic case, the 
number of bridges m is a Poisson distributed random 



number P m = rn" 1 exp (— m) /ml where m denotes the 
mean of the distribution. After staying for some time at 
a particular site, the particle either breaks all its bridges 
and departs, or hops a distance a to a new site charac- 
terized by a new value for the number of bridges m. In 
this sense we have coarse-grained the particle motion by 
discrete steps of the correlation length a, the distance 
after which the number of bridges becomes statistically 
independent of the value at the previous location. 




FIG. 1: (Color online). A snapshot of particles interacting 
with a two-dimensional substrate. Particles are alternately 
bound to the substrate by bridges, or unbound and free to 
diffuse in a direction normal to the substrate plane. 



There is a percolation transition which separates the 
diffusive regime from the localized regime. In the lo- 
calized regime the particle remains close to the original 
location until breaking all its bridges and departing. In 
the diffusive regime the particle undergoes a random walk 
by breaking and reforming multiple bridges. The tran- 
sition between the two regimes occurs at the percolation 
threshold where one first encounters an infinitely con- 
nected cluster of sites with m > 0. Since the critical 
probability of bond percolation on the square lattice (l3| 
is i, the transition is given by Pq — \. Below, we cal- 



2 



culate the departure time distribution $(i) in both the 
localized and diffusive regimes, and study the random 
walk statistics above the percolation threshold. 

Localized regime. Consider a particle attached 
to the substrate by to bridges below the percolation 
threshold m = log 2. The probability that the parti- 
cle departs from the surface between time t and t + dt 
is $ m (t)dt ~ K m exp [— K m t] dt. The departure rate 
K m — ^ exp (—Am) is given by the Arrhenius relation 

K m ~ exp ^ fcjjr ) with tq a characteristic timescale for 
bridge formation. By averaging this distribution over 
the statistics of m we arrive at: 



(1) 



m— 1 



When performing the averaging we do not include states 
with m = bridges. For this reason we work with 
a renormalized probability distribution P m — P m /(1 — 



exp (— m)) so that 



P — 1 
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FIG. 2: (Color online). Departure time distribution function 
vs. time at the percolation threshold m = log 2. The solid 
lines are calculated from Eq. [T]in the localized regime, and 
the dotted curves are calculated from Eq. [4] in the diffusive 
regime. 

The results of this calculation can be compared to a re- 
cent experiment which determined the time-varying sep- 
aration of two DNA-grafted particles in an optical trap 
[l[ . In this experiment two particles are bound by DNA 
bridges, and after breaking the connections diffuse to the 
width of the optical trap. Because the length of the DNA 
chains grafted on the particle is much less than the par- 
ticle radius, surface curvature effects can be neglected. 
Hence, the experiment resembles a particle interacting 
with a localized site on a 2D substrate. Experimentally, 
the tail of the departure time distribution was observed 
to be a power law <I>(i) ~ t^ 1 - 5 . Such behavior is in- 
deed reproduced by Eq. [T]with a binding free energy on 



the order of several fcgT. For strong enough binding, 
A > 1, the departure time distribution function exhibits 
a pronounced multi-exponential character. 

Diffusive regime. The departure time distribution 
is strongly altered above the percolation threshold. In 
this regime the particle can move around to find a more 
favorable state on the surface. This leads to a much 
longer lifetime of the bound state, a phenomenon sim- 
ilar to aging in glassy systems. The hopping rate be- 
tween two neighboring sites is given by the Arrhenius 
law, /tj_>j ~ exp[—A(m,i — rrij)8 (rrii — rrij)]. Here 
0(x) is the Heaviside step function. The problem can be 
greatly simplified since the ensemble averaged hopping 
rate from a site with to bridges can be well approximated 
by an effective Arrhenius relation: 



K m = — exp [— A(m — m)] 
to 



(2) 



In the case when Am is sufficiently large, the 

probability of staying attached to the surface af- 

/ „ \n-l 

ter an n step random walk is 1 1 — 1 = 

[1 — exp (— Am)]" 1 . Interestingly, this probability is 
virtually independent of the particular bridge numbers 
{mi,...,m„} realized during the walk. We conclude 
that the probability of particle departure after exactly 
n steps is /„ = [exp (7) — I] exp(— 771), where 7 = 
— log [1 — exp (— Am)]. If we let (j) n (t) denote the de- 
parture time distribution for a walk of n steps we have: 
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To complete the calculation it is convenient to Fourier 
transform <p n (t). One can then sum the resulting ge- 
ometric series to obtain $>(lu) and perform the inverse 
transform to derive the following result: 

= exp( 7 ) [exp( 7 ) - 1] £ ex ^~ z f (4) 
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Here z r labels the roots of the equation 
exp 7 ~ >. 



0. 
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In Figure [5] the two results for the departure time 
distribution are compared at the percolation threshold 
to = log 2. For fixed to a change in A is directly re- 
lated to a change in the average binding free energy. As 



FIG. 3: (Color online). Departure time distribution function 
vs. time as determined by Eq. [4] in the diffusive regime. In 
the plot Am = 4. 



FIG. 4: (Color online). Rms displacement vs. time holding 
Am = 4. The inset is the dimensionless diffusion coefficient 
defined as D = plotted against time. 



indicated in the plot, increasing A decreases the rate of 
particle departure. 

In Figure [3] we plot the departure time distribution as 
determined by Eq. |4] in the diffusive regime. In the 
figure we hold the product Am = 4 constant. The op- 
timal regime for fast departure is to have a large num- 
ber of weakly bound bridges. In this case the depar- 
ture time is well approximated as a single exponential, 
$(*) = K w ex^{-K w t). 

Finally, we discuss the statistics of the in-plane diffu- 
sion of the particle. We notice that the in-plane trajec- 
tory of the particle subjected to a delta-correlated ran- 
dom potential remains statistically equivalent to an un- 
biased random work. Therefore, the mean-squared dis- 
placement after n steps is given by (r 2 ^ = na 2 . However, 
as the particle explores the landscape the average hop- 
ping time becomes longer and the diffusion gets slower. 
In the limit n — > oo, the average hopping time can be 
determined from the equilibrium canonical distribution. 
For the case of Poisson distributed m, this corresponds 
to a finite yet renormalized diffusion coefficient D* with 
D = a 2 /4r : 
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exp (Am) [exp (m) — 1] 
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However, it may take a very long time to achieve this 
"ergodic" behavior. In the transient regime, an n— step 
random walk cannot typically visit sites with an arbi- 
trarily large number of bridges m. Instead, one should 
average the hopping times only over sites with m < m*. 
In the language of the statistics of extreme events, m* — 1 
is the maximum " expected" value of m in a sample of n 



independent events 



scpec 



Both the average diffusion time 



expressed in terms of m* , which defines their relationship 
in parametric form: 



(r 2 ) = - 



(t) 



(r 2 ) 



D* 



P(m, m*) 



P(me j m*) 
1 — exp (— me A ) 
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(i), and mean square displacement (r 2 ^ = na 2 , can be 



Here P(x,m*) = j(x, m*)/T (m*) = 

oo 

exp (— x) ^2 x k /k\ is the regularized lower incom- 

k— m* 

plete r function. It is easy to see that in the limit 
to* — > oo we recover the renormalized diffusion relation 
(t) = (r 2 ) /D*, although this occurs at very long, often 
unrealistic times. In the transient regime we expect 
anomalous, subdiffusive behavior. As shown in Figure 
IU this regime is typical for strong enough key-lock 
interactions. The predicted anomalous diffusion may be 
well described by a power law, (r 2 (i)) ~ (t) a , with a 
non-universal exponent a < 1 . 

This work provides additional insight into the slow 
crystallization dynamics of key-locking particles (see fig- 
ure 1 /im diameter particles were grafted 
with ssDNA and formed reversible, disordered aggre- 
gates. The average number of key-lock bridges between 
particle pairs was TO ~ 1 . By further reducing the graft- 
ing density of DNA strands on the particles (m < 1), the 
authors of reference [J] observed random hexagonal close- 
packed crystals. Crystallization requires that colloids 
repeatedly depart and reattach to the growing structure, 
in an effort to find their desired lattice location. We can 
attempt to quantify this optimal experimental regime of 
fast departure by determining the time T required for 
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FIG. 5: (Color online). Schematic depiction of key-lock bind- 
ing between nanoparticles functionalized with complementary 
ssDNA. The resulting structures can be disordered, fractal- 
like aggregates, or crystalline. 



90% of the particles to depart: 0.1 = $>(t)dt. Fig- 



ure [6 

ergy 



is a plot of T vs. m at constant binding free en- 
+ log(l — exp(— m)) = const. The 
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optimal regime is to have a large number (m ~ 10) of 
weakly bound bridges. To realize this regime exper- 
imentally we propose the introduction of long, flexible 
DNA linkers to a system of particles with a high cover- 
age of short ssDNA. This scheme increases the number 
of key-lock bridges between particle pairs as compared to 
previous experiments, and therefore has the potential to 
substantially reduce the time required for crystallization. 
The prediction of a localized regime below the percolation 
threshold (the local minima in Figure [6]) where particle 
departure is relatively fast is confirmed experimentally 
by the crystallization observed in [l[. 




FIG. 6: (Color online). Plot of the time T required for 90% 
of the particles to depart vs. m. Note the jump in T which 
occurs at the percolation threshold separating the localized 
regime from the diffusive regime. 

In this work we studied the dynamics of particles which 



form multiple, reversible key-lock bridges. There is 
a percolation transition which separates the regime in 
which particles are localized near their original location 
from the regime where they exhibit diffusive behavior 
by breaking and reforming bridges. At low coverage 
the key-locking system exhibits a broad, power law like 
distribution of departure times, but no lateral diffusion. 
Above the percolation transition (m = log 2) diffusion 
allows the particle to cascade into deeper energy wells 
with a large number of key-lock bridges. This leads to 
an increase in the bound state lifetime similar to aging in 
glassy systems. The statistics for the particles' in-plane 
diffusion were determined. For relatively weak key-lock 
interactions there is a finite renormalization of the dif- 
fusion coefficient. However, as A increases, the system 
exhibits anomalous, sub-diffusive behavior analogous to 
dispersive transport in disordered semiconductors. We 
discussed the connection between our work and recent ex- 
periments with DNA-coated colloids. The findings indi- 
cate that the optimal regime for colloidal crystallization, 
where particle departure is a relatively fast process, is to 
have a large number of weakly bound key-lock bridges. 
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